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Abstract 

Beginning with the general Bjorken-Mtingwa solution for intrabeam scattering (IBS) we derive 
an accurate, greatly simplified model of IBS, valid for high energy beams in normal storage ring 
lattices. In addition, we show that, under the same conditions, a modified version of Piwinski's 
IBS formulation (where r]1y/ I3x,y has been replaced by 'Hx,y) asymptotically approaches the result 
of Bjorken-Mtingwa. 
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INTRODUCTION 



Intrabeam scattering (IBS), an effect that tends to increase the beam emittance, is impor- 
tant in hadronic|l[] and heavy ion[3 circular machines, as well as in low emittance electron 
storage rings]^. In the former type of machines it results in emittances that continually 
increase with time; in the latter type, in steady-state emittances that are larger than those 
given by quantum excitation/synchrotron radiation alone. 

The theory of intrabeam scattering for accelerators was first developed by PiwinskiQ, a 
result that was extended by MartinijQ, to give a formulation that we call here the standard 
Piwinski (P) methodj^; this was followed by the equally detailed Bjorken and Mtingwa 
(B-M) result 0. Both approaches solve the local, two-particle Coulomb scattering problem 
for (six-dimensional) Gaussian, uncoupled beams, but the two results appear to be different; 
of the two, the B-M result is thought to be the more general]^. 

For both the P and the B-M methods solving for the IBS growth rates is time con- 
suming, involving, at each time (or iteration) step, a numerical integration at every lattice 
element. Therefore, simpler, more approximate formulations of IBS have been developed 
over the years: there are approximate solutions of ParzenQ, Le Duff|Ty], Raubenheimer |IT| , 



and Wei||12||. In the present report we derive — starting with the general B-M formalism — 
another approximation, one accurate and valid for high energy beams in normal storage ring 
lattices. We, in addition, demonstrate that under these same conditions a modified version 
of Piwinski's IBS formulation asymptotically becomes equal to this result. 

HIGH ENERGY APPROXIMATION TO BJORKEN-MTINGWA 
The General B-M Solution[0 

Let us consider first machines with bunched beams that are uncoupled and have vertical 
dispersion due to e.g. orbit errors. Let the intrabeam scattering growth rates be defined as 
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with (Tp the relative energy spread, ex the horizontal emittance, and ey the vertical emittance. 
The growth rates according to Bjorken-Mtingwa (including a -\/2 correction factor |jl3|, and 
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including vertical dispersion) are 
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where i represents p, x, or y; 



[det(L + A/)]V2 
1 \ u.f 1 



L + \I 



A 



rlcN 



L + XI 



(2) 



(3) 



with tq = 2.82 X 10^^^ m, the classical electron radius, c the speed of light, the bunch 
population, /3 the velocity over c, 7 the Lorentz energy factor, and as the bunch length; (log) 
represents the Coulomb log factor, () means that the enclosed quantities, combinations of 
beam parameters and lattice properties, are averaged around the entire ring; det and Tr 
signify, respectively, the determinant and the trace of a matrix, and / is the unit matrix. 
Auxiliary matrices are defined as 
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The dispersion invariant is 7i = [77^ + (Prj' — \[3'rjY]/ [3, and = 77' — ^P'r]/P, where (3 and 
7] are the beta and dispersion lattice functions. 

For unbunched beams ag in Eq. || is replaced by C/{2\/2n), with C the circumference of 
the machine. 
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The Bjorken-Mtingwa Solution at High Energies 



Let us first consider 1/Tp as given by Eq. |^. We first notice that, for normal storage ring 
lattices (where {'Hx,y/Px,y) ^ 1), the off-diagonal elements in L, —70, are small and can be 
set to zero. Then all matrices are diagonal. Let us also limit consideration to high energies, 
i.e. let us assume a,b -C 1, with 
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with 

^ = ^ + ^ + ^ . (9) 

O'h CTp ^x £y 

Note that if a,b <^ 1, then the beam is cooler longitudinally than transversely. If we consider, 
for example, KEK's ATF, a 1.4 GeV, low emittance electron damping ring, ey/ex ~ 0.01, 
a -0.01, 6 -0.10. 

If the high energy conditions are met then the 2nd term in the braces of Eq. ^ is small 
compared to the first term, and can be dropped. Now note that L2,2 can be written as '-f'^/ajj. 



For high energy beams a factor in the denominator of the integrand of Eq. ^ yl'^/'^H + -^j 
can be approximated by 'y/crn', also, the (2,2) contribution to Tr[(L + A/)^^] becomes small, 
and can be set to 0. Finally, the first of Eqs. |] becomes 
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A plot of g{a) over the interval [0 < a < 1] is given in Fig. |l|; to obtain the results for a > 1, 
note that g{a) = g{l/a). A fit to g, 

(7(a) ^ 2a(°-°2i-°-°^^i'^°) [for O.OK a < 1] , (12) 

is given by the dashes in Fig. |T|. The fit has a maximum error of 1.5% over [0.02 < a < 1]. 
Similarly, beginning with the 2nd and 3rd of Eqs. ^, we obtain 
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FIG. 1: The auxiliary function g{a) (solid curve) and an analytical approximation, g 

2^,(0.021-0.044 Ina) (dashes). 



Our approximate IBS solution is Eqs. p!0|J13| . Note that Parzen's high energy formula is a 
similar, though more approximate, result to that given hereQ; and Raubenheimer's approx- 
imation is formulas similar, though less accurate, than Eq. 10 and identical to Eqs. 1^[11 



Note that the beam properties in Eqs. |T0|,|T3|, need to be the self-consistent values. Thus, 
for example, to find the steady-state growth rates in electron machines, iteration will be 
required. Note also that these equations assume that the zero-current vertical emittance is 
due mainly to vertical dispersion caused by orbit errors; if it is due mainly to (weak) x-y 
coupling we let Hy = 0, drop the l/Ty equation, and simply let ey = ne^, with k the coupling 
fact or 101 . 



COMPARISON TO THE PIWINSKI SOLUTION 



The Standard Piwinski Solution H 



The standard Piwinski solution is 
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Parameters are: 
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The function / is given by: 
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where 



(18) 



The parameter d functions as a maximum impact parameter, and is normally taken as the 
vertical beam size. 

Comparison of Modified Piwinski to the B-M Solution at High Energies 

To compare with the B-M solution, let us consider a slightly changed version of Piwinski 
that we call the modified Piwinski solution. It is the standard version of Piwinski, but with 
rj'^/P replaced by Ti. {i.e. d, b, ah, become a, b, an, respectively). Let us also assume high 
energy beams, i.e. let a,b <^ 1. 

Let us sketch the derivation. First, notice that in the integral of the auxiliary function / 
(Eq. |l^): the —0.577 can be replaced by 0; the —3u^ in the numerator can be set to 0; P (Q) 
can be replaced by Va^ + (y/W^-v^). The first term in the braces can be approximated 
by a constant and then be pulled out of the integral; it becomes the effective Coulomb log 
factor. Note that for the proper choice of the Piwinski parameter d, the effective Coulomb 
log can be made the same as the B-M parameter (log). For flat beams (a <^ 6), the Coulomb 
log of Piwinski becomes (log) = In [da'jj/{Aroa'^)]. 

We finally obtain 
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The integral is an elliptic integral. The first of Eqs. then becomes 
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with 
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We see that the the approximate equation for 1/Tp for high energy beams according to 
modified Piwinski is the same as that for B-M, except that h{a,b) replaces g{a/b). 

We can now show that, for high energy beams, h{a,b) ^ g{a/b): Consider the function 
h{a, b, (), which is the same as h{a, b) except that the upper hmit of integration is infinity, and 
the in the denominator are replaced by (u"^. It is simple to show that d(^h{a,b,()\c=i = 
g{a/b) = h{a,b,l). Now for high energies {a,b small), reducing the upper limit in the 
integral of h{a,b, 1) to 1 does not significantly change the result, and h{a,b) ~ g{a/b). To 
demonstrate this, we plot in Fig. ^ the ratio h{a,b)/g{a/b) for several values of a. We 
see, for example, for the ATF with ey/ex ~ 0.01, a ~ 0.01, a/b ~ 0.1, and therefore 
h{a,b)/ g{a/b) = 0.97; the agreement is quite good. 
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FIG. 2: The ratio h[a,h) / g[a/h) as function of a/6, for three values of a. 



Finally, for the relation between the transverse to longitudinal growth rates according to 
modified Piwinski: note that for non-zero vertical dispersion the second term in the brackets 
of Eqs. |14| (but with rjl y/ Px,y replaced by Hx^y) , will tend to dominate over the first term, 
and the results become the same as for the B-M method. 

In summary, we have shown that for high energy beams {a,b -C 1), in rings with a 
standard type of storage ring lattice: if the parameter d in P is chosen to give the same 
equivalent Coulomb log as in B-M, then the modified Piwinski solution agrees with the 
Bjorken-Mtingwa solution. 
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NUMERICAL COMPARISON||3| 



We consider a numerical comparison between results of the general B-M method, the 



modified Piwinski method, and Eqs. p!0|JT3| . The example is the ATF ring with no coupling 
and vertical dispersion due to random orbit errors. For our example = 17 fim, yielding 
a zero-current emittance ratio of 0.7%; the beam current is 3.1 mA. The steady-state growth 
rates according to the 3 methods are given in Table 1. We note that the Piwinski results are 
4.5% low, and the results of Eqs. rU,0|, agree very well with those of B-M. Finally note that, 
not only the growth rates, but even the differential growth rates — i.e. the growth rates as 
function of position along the ring — agree well for the three cases. 



TABLE L Steady-state IBS growth rates for an ATF example including vertical dispersion due to 
random errors. 



Method 



1/T, [s-i] 



[s-i] 



1/Ty [s-i] 



Modified Piwinski 
Bj orken-Mtingwa 



Eqs. iq,l 



25.9 
27.0 
27.4 



24.7 
26.0 
26.0 



18.5 
19.4 
19.4 
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